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Charge order and superconductivity in a two-dimensional triangular lattice at = 2/3 
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To investigate the possibility of charge order and superconductivity in a doped two- 
dimensional triangular lattice, we study the extended Hubbard model with variational Monte 
Carlo method. At n = 2/3, a commensurate filling for a triangular lattice, it is shown that the 
nearest-neighbor Coulomb interaction V induces honeycomb-type charge order and antiferro- 
magnetic spin order at U > lOt. We also discuss the possibility of superconductivity induced 
by charge fluctuation and the relation to the superconductivity in Nao.35Co02 ■ I.3H2O and 
6'-type organic condoctors. 
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Since the discovery of high- Tc cuprates, it have 
been clarified that superconductivity and magnetism are 
closely connected with each other in strongly correlated 
electron systems. Hubbard model and t-J model in a two- 
dimensional square lattice have successfully described 
the relation between them and showed that antiferro- 
magnetic spin fluctuation induces dj.2_j^2-wave supercon- 
ductivity in doped Mott insulators. 

In a frustrated system like a two-dimensional triangu- 
lar lattice, on the other hand, the relation between su- 
perconductivity and magnetism have not been clarified 
yet. Generally, frustration suppresses long-range mag- 
netic orders and induces large fluctuation. Anderson flrst 
proposed a resonating valence bond (RVB) state for the 
ground state of an antiferromagnetic two-dimensional tri- 
angular lattice, as a notion substituted for the Fermi liq- 
uid. ""^ Since then, study of frustrated spin systems has 
been greatly improved and various interesting phenom- 
ena have been clarified. However, when some carriers are 
doped into the frustrated lattices, the character of elec- 
tron systems has not been clarified extensively. Actually, 
materials considered as doped frustrated systems like 9- 
type organic conductors^ and Naa;Co02 • yH20^''' show 
various magnetic and charge orders, as well as unconven- 
tional superconductivity. 

In a system away from half filling, carriers obtain a 
large mobility and the charge fluctuation is enhanced by 
strong electron correlation. Then, the charge degrees of 
freedom become important as well as the spin degrees of 
freedom. To deal with the charge fluctuation correctly, 
we have to take into account the long-range Coulomb in- 
teraction which is not contained in the Hubbard model 
or t-J model. An extended Hubbard model with on-site 
Coulomb interaction U and nearest-neighbor one V is 
an approriate model to treat this problem. The compe- 
tition between U and V is considered to induce vari- 
ous states. For example, in a two-dimensional square lat- 
tice with quarter-filled band, previous theoretical studies 
show that a checkerboard-type charge order (CO) be- 
comes stable at large V region^ and rf^jy-wave supercon- 
ductivity is stabilized at intermediate V region.^' ^ 




Fig. 1. First Brillouin zone and Fermi surface used in this paper. 

In this letter, we study an extended Hubbard model 
in a two-dimensional triangular lattice with variational 
Monte Carlo method having the Na:rCo02 system in 
mind. Although the multi-orbital effect may be impor- 
tant for superconductivity,^'^ we focus on the single- 
band model^°"^^ with charge fluctuation. ^^"^^ 

We show that the competition between U and V in- 
duces the two types of CO states and discuss the possibil- 
ity of superconductivity induced by charge fluctuation. 

We introduce the following extended Hubbard model 
in a two-dimensional triangular lattice, 

-ff = ^ (^kcl^Cka + U'^ni^nii-\-V'^ n^rij , (1) 

fccr i 

where denotes the summation over the nearest- 

neighbor sites. The bare energy dispersion is given by 

1 

Efc = — 2t(cos fca; -t- 2 cos -fcj, cos -^fcy). (2) 

In this non-bipartite lattice, there is no electron-hole 
symmetry and the sign of t is crucial. The sign of t should 
be determined so as to reproduces the shape of the Fermi 
surface obtained from experiments or band calculations 
for each material under consideration. For the supercon- 
ducting material Nao.35CoO2-l.3H2O, the choice oit <0 
and n = 4/3 reproduces the large Fermi surface around F 
point calculated by the LDA calculation^^ or observed by 
ARPES experiments.^"' In the following, we use t > 
and n = 2/3 by using the electron- hole transformation 
as shown in Fig. 1. To study the ground state, we use 
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Fig. 2. Three different cliarge and spin ordered states in a two- 
dimensional triangular lattice for n = 2/3. (a) Paramagnetic CO 
where solid (open) circles represent preferred B and C (avoided 
A) sublatticcs. (b) Bipolaronic state where electrons prefer A 
sublattice in spite of double occupancy, (c) AF - Antiferromag- 
netic state with honeycomb-type charge order, (d) The folded 
first Brillouin zone for a honeycomb-type charge order. Dashed 
line represents the original first Brillouin zone. 



variational Monte Carlo (VMC) method. We introduce 
the Jastrow-type trial wave function for the extended 
Hubbard model, 



where 



I*) = pwPg m , 

i 



(3) 

(4) 

(5) 



Pg is a Gutzwiller projection operator which reduces the 
probability of double occupancy of electrons at the same 
site. Pw is a nearest-neighbor projection operator which 
controls the weight of the electron configuration at the 
nearest-neighbor sites. This kind of projection operator 
is successfully used in the square lattice. ^^'^"^ g and w 
are variational parameters with < ^f, tz; < 1. If we set 
gf = w = 1, l^*) is reduced to |<I>) which is constructed 
from a mean-field solution of a certain Hamiltonian and 
its explicit form is discussed shortly. 

In a two-dimensional triangular lattice at n = 2/3, 
several honeycomb-type states as shown in Fig. 2 are 
commensurate charge-ordered states. The competition 
between U and V is expected to induce some charge- 
ordered states. In the following, we denote the three dif- 
ferent states as paramagnetic charge order (PCO), bipo- 
laronic charge order (BP), and antiferromagnetic charge 
order (AF) as shown in Fig. 2(a)-(c). Trial wave functions 
corresponding to these states are obtained from the so- 



lution of the mean-fileld Hamiltonian 

-ffa = X] ^kci^Cka + 2Ac ^ cos(Q • ri)ni 



ka 



2As ^ cos (Q-ri + 



(6) 



where Ac is an order parameter of charge order, Ag is an 
order parameter of spin order and Q = 0) is an or- 
dering wave vector. PCO is a state in which electrons pre- 
fer B and C sublattices and avoid A sublattice, and there 
is no spin order. Thus we can choose Ac > 0, As = 0. BP 
is a bipolaronic state in which electrons prefer A sub- 
lattice and avoid B and C sublattices with double occu- 
pancy, i.e., Ac < 0, As = 0. Finally, AF is a state that 
charge order and antiferromagnetic spin order coexist^'' 
(Ac > 0, As > 0). If we set Ac = Ag = 0, the state 
reduces to paramagnetic metal (PM). 

For the trial wave fuction of superconductivity, we em- 
ploy the BCS type wave function, 

\^)=J{{uk+Vkci^cl^^)\Q). (7) 

k 

The form of this wave function is not convenient for the 
VMC method because the number of electrons are not 
fixed. To obtain the wave function with a fixed number 
of electrons, we extract the term imvolving A^e electrons 

as,^^ 



iVe/2 

|$)cx I ^a,,cl^c]J |0) 



(8) 



where 



1 



E 

k 

E- 



^gifc-(rj-rj) 
Uk 



Afc 



^ik-{r 



(9) 



The matrix element refiects the symmetry of Cooper 
pairs: for singlet pairing, a 



ttij and for triplet. 

The possible paring symmetry of A^ , the gap 
function of superconductivity, is determined by group 
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Fig. 4. Band dispersion obtained from the mean-field solution of 
eq.(6). (a) Paramagnetic metal (Ac = 0), (b) charge ordered 
metal (Ac = — O.lt), (c) charge ordered insulator (Ac = —0.5t) 
which is the BPI. Horizontal lines represent the Fermi energy. 

theory. The six candidates are s,d,i (singlet) and p,f, 
next-nearest-neighbor / (triplet) symmetries. 

We determine the phase diagram in a,U — V plane by 
optimizing the ground state energies of trial wave func- 
tions. Figure 3 shows the phase diagram in a U—V plane. 
For large value of V, BPI ("I" denotes insulator) state 
is stable because the effect of V overcomes U and the 
electrons concentrate in one of the three sublatticcs in 
spite of double occupancies. At U > lOt, there appears 
an AFT region between paramagnetic metal (PM) and 
BPI regions. For small U, the state changes as PM 
BPI as V increases. For large U, however, electrons try 
to order without loss of on-site potential energy. Then 
the state changes as PM AFI BPI. Since the en- 
ergy loss in the Coulomb interactions in BPI and AFI 
arc f7/3 and V per site, the critical value of V ~ U/3 
can bo easily expected. On the other hand, we find that 
PCO region does not exist in this phase diagram because 
the calculated total energy of PCO is always higher than 
AFI. That is to say, the order parameters of charge order 
and spin order grow simultaneously. 

Interestingly, we find that all three phase transitions 
(PM BPI, PM AFI and AFI BPI) are first order. 
At V = Vc, the critical value of V, the order parameter 
has a jump and the variational energies of the states 
cross. Reflecting the first order character, the physical 
quantities like magnetization or charge density change 
discontinuously. In order to understand the reason of this 
first order transition, we show the change of the band 
structure in Fig. 4. When Ac = 0, the system is a para- 
magnetic metal (PM, Fig. 4(a)) state. As Ac grows, two 
degenerate bands along A-B line (see Fig. 2(d)) begin 
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Fig. 5. Schematic representation of next-nearest-neighbor /-wave 
symmetry in a fc-space and a real space. A^ changes its sign 
six times in a fc-space and it corresponds to the next-nearest- 
neighbor paring in a real space. 

to split. When |Ac| < t/3, the system shows a metallic 
behavior because there is an overlap between the upper 
band and lower band. We call this state as "charge or- 
dered metal (COM, Fig. 4(b))" . If |Ac| > t/3, the overlap 
disappears and we call this state as " charge ordered in- 
sulator (CGI, Fig. 4(c))". PM BPI transiton in Fig. 3 
is, that is to say, "PM to COI" type first-order transi- 
tion. COM state is not realized in this condition. This 
result is interpreted as follows. Since n = 2/3 is a com- 
mensurate filling for a two-dimensional triangular lat- 
tice, the honeycomb-type charge orders like BPI or AFI 
are greatly favored. We consider that the optimal carrier 
doping removes the frustration and reduces the charge 
fluctuation. COM is an intermediate state with weak sta- 
bility and thus it is suppressed in this condition. There- 
fore, if we slightly break the commensurability, for exam- 
ple, by the introduction of a lattice anisotropy or doping, 
the fluctuation is enhanced again and COM state can 
appear with second-order transition. Indeed, the study 
of the extended Hubbard model in a two-dimensional 
anisotropic triangular lattice by Seo shows the existence 
of COM state in some parameter region. 

Finally, let us discuss the possibility of superconduc- 
tivity in the present case. Our calculation shows that 
there is no region where a superconductivity becomes 
stable. This result is reasonable at least for = be- 
cause the nesting of Fermi surface is absent. Further- 
more, the van Hove singularity which is advantageous 
to superconductivity is away from the Fermi surface. 
The third-order perturbation calculation^'' and the one- 
roop renormalization-group calculation^* for the Hub- 
bard model show the same results. Then the problem is 
whether the introduction of V induces a superconductiv- 
ity or not. Qualitatively speaking, possible candidate is 
next-nearest-neighbor (nnn) /-wave pairing. When V is 
large, electrons avoid each other and nnn pairing is con- 
sidered to be favored in the real-space representation. 
Schematic explanation is shown in Fig. 5. Indeed, we 
flnd that nnn /-wave pairing has a lower energy than the 
paramagnetic metal state at large value of V. However, 
its condensation energy is quite small compared with 
those of charge ordered states (BPI and AFI). Motrunich 
and Lee^^ also showed that the obtained is much 
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smaller than the realistic value, although the most stable 
pairing is nnn /-wave near the charge order instability. 
We specurate that the reason why the superconductivity 
is not stabilized is the same as the case of COM. Thus, 
the appropriate breaking of the commensurability will 
induce the large fluctuation and then the charge order 
and superconductivity can compete each other. 

In the following, we propose two ways to break the 
commensurability and induce the superconductivity. One 
is to change the electron density away from the commen- 
surate filling. The RPA study of the extended Hubbard 
model in a two-dimensional triangular lattice at n = 0.8 
and 1.2 by Tanaka et al}^ shows that nnn /-wave su- 
perconductivity is realized near the CDW instability. In 
this condition, CDW state is incommensurate and its 
stability is not so strong. Therefore, the charge fluctu- 
ation is enhanced near this instability and nnn /-wave 
superconductivity appears. Since the electron filling of 
Nao.35Co02 • I.3H2O has not been determined precisely 
because of the existence of H3O+ ion, the above mecha- 
nism might be applicable to this material. 

Another way is to introduce some anisotropy. Super- 
conductivity in ^-type organic conductors are located 
next to the charge order state and its origin is consid- 
ered to be the charge fluctuation. They have quasi-two- 
dimensional anisotropic triangular lattice with n = 1.5 
filling. For the square lattice, this filling is commensu- 
rate and leads to the checkerboard-type charge order. In 
real materials, however, the anisotropics of hopping or 
off-site Coulomb interaction break this commensurabil- 
ity and induces various states like several charge orders 
(three kind of stripe statcs^^ and 3-fold state^°) and, 
probably, unconventional superconductivity. We can re- 
late this mechanism to the one discussed so far. The 
study of an anisotropic triangular lattice will interpo- 
late the square- and isotropic triangular lattice and give 
the consensus about the charge fluctuation. It is left for 
the future problem. 

In summary, we have studied the extended Hubbard 
model in a two-dimensional triangular lattice and find 
the two types of charge order states (BPI and AFI) in 
a, U — V phase diagram. These states are so stable at 
n = 2/3 that the superconductivity induced by charge 
fiuctuation has no chance to appear. However, the ap- 
propriate breaking of the commensurability can induce 
the charge fluctuation and the possibility of the super- 
conductivity is expected to arise. We expect that this 
mechanism is applicable to ^-type organic superconduc- 
tors or Nao.35Co02 • I.3H2O. 
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